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Abstract 


Phase sequences of the biaxial nematic liquid crystal in the interior 
of the essential triangle are studied with Wang Landau sampling. The 
evidence points to the existence of an intermediate unixial phase with 
low biaxiality in the isotropic to biaxial nematic phase sequence. 
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1 Introduction 


The biaxial liquid crystal phase (iV^lwhich was predicted in the early 70’s 
[1] continues to be elusive inspite of signihcant progress made in the theoretical [2] 
- [E], experimental [T5] - [20] and computer simulation [2T] - [27] studies re¬ 
garding the existence and properties of the phase. The recent mean held the¬ 
oretical (MFT) studies of a quadrupolar Hamiltonian model [5] - [9] predict 
a universal mean held phase diagram for biaxial nematic along the boundary 


of a triangular parameter space OIV( see Fig, 1(a)) wherein the condensa¬ 


tion of the biaxial phase could occur either from the uniaxial (Nu) phase or 
directly from the isotropic phase (/). These predictions, which were partly 
verihed by Monte Carlo simulations, were found to be unsatisfactory in the 
limit of vanishing biaxial-biaxial interaction in the repulsive region for the 
Hamiltonian, thus requiring further study. 

In this context, our recent WL simulations of the phase sequences along 
the boundary of the triangle OIV [2Sl 121] suggested a qualitative mod- 
ihcation of the MFT phase diagram as the Hamiltonian is driven to the 
partly repulsive regions. The efficient entropic sampling technique employed 
seeking otherwise inaccessible rare microstates, pointed to the 
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Figure 1: (a)The essential triangle, depicting the interior trajectory IZW 
along which detailed simulations have been carried out. (b)Location of des¬ 
ignated values of ( 7 , A) along IW (also see Table. ([I])) 


existence of possible hindering free energy barriers within the system result¬ 
ing from the absence of stabilising long range order of one of the molecular 
axes. Keeping in view the crucial role played by the degree of cross-coupling 
between the uniaxial and biaxial tensorial components of the neighbouring 
molecules in the condensation of the biaxial phase, we present in this pa¬ 
per, the results of a similar detailed simulation study which was carried out 
along a segment IW in the interior of the essential triangle, where W is the 
midpoint of OV (see Fig. [1]). 

We carried out a systematic simulation study using the entropic sampling 
technique (WL algorithm) along the segment IW in order to obtain a generic 
phase diagram inside the essential triangle. The mean field Hamiltonian and 
the simulation model are presented in section II. The sampling technique and 
details of simulation are discussed in section III. The observations from the 
simulations are reported in section IV and the conclusions are discussed in 
section V. 


2 Hamiltonian model 

The MF analysis [S] - [H], is based on the general quadrupolar orientational 
Hamiltonian, proposed by Straley [2] and set in terms of tensors [S]. Ac¬ 
cordingly, the interacting biaxial molecules are represented by two pairs of 
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symmetric, traceless tensors (q, b) and (q\ b'). Here q and q' are uniaxial 
components about the unit molecular vectors m and m , whereas b and b 
(orthogonal to q and q , respectively), are biaxial. These irreducible compo¬ 
nents of the anisotropic parts of susceptibility tensor are represented in its 
eigen frame (e, ej_, m) as 


q ■= m® m - (la) 

3 

6 := e (g) e — ej_ (g) ej_ (lb) 

where I is the identity tensor. Similar representations hold for q and b' in 
the eigen frame (e , e_|_, m ). The interaction energy is written as 


H = -U[i qq' + -i{qb' + q' ■b) + \bb'] (2) 

where U is the scale of energy, ^ = ± 1 , 7 and A are dimensionless interac¬ 
tion parameters, determining the relative importance of the uniaxial-biaxial 
coupling and biaxial-biaxial coupling interactions between the molecules, re¬ 
spectively. 

Mean-held analysis of the Hamiltonian identihes a triangular region OIV 
in the ( 7 , A) plane - called the essential triangle - representing the domain 
of stability into which any physical system represented by Eqn. ([2]) can be 
mapped mu (see Fig. [1]). The line C 1 C 3 is a tricritical line whereas C 2 C^ is a 
triple line. The dispersion parabola A = traverses through the interior 

of the triangle, intersecting IV at the point T, called the Landau point. 
Region of the triangle above the parabola corresponds to a Hamiltonian 
where all the terms are attractive, while the region below is noted to be 
partly repulsive [7]. MFT predicts a / —>■ Nu Nb phase sequence in the 
quadrangle OC 2 C 3 V and a direct I —)• Nb transition in the triangular region 
IC 2 C,. 

For simulation purposes, the general Hamiltonian in Eqn. ([2]) is conve¬ 
niently recast as a biaxial mesogenic lattice model, where particles of D 2 h 
symmetry, represented by unit vectors Ua, Vb on lattice sites a and b interact 
through a nearest-neighbour pair potential [33] 


U = -e{Gs 3 - 27 (Gii - G22) + A[ 2 (Gn + G22) - G33]}. ( 3 ) 

Here fab= (Va-Ub), Gab=P 2 {fab) with P 2 denoting the second Legendre 
polynomial. The constant e (set to unity in simulations) is a positive quan- 
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tity setting the reduced temperature = ksT/e, where T is the absolute 
temperature of the system. 

3 Details of Simulation 

The Wang-Landau (WL) sampling [30] is a flat histogram technique designed 
to overcome energy barriers encountered, for example, near first order tran¬ 
sitions, by facilitating a uniform random walk along the energy {E) axis 
through an appropriate algorithmic guidance. The sampling, originally de¬ 
veloped for Hamiltonian models involving random walks in discrete conhgura- 
tional space, continues to be applied to various problems in statistical physics 
[5n 133] , polymer and protein studies |3S1 ETl [3S] and is being developed for 
more robust applications for continuous systems [3H1 SDl IHl 1121 SSI SH ES] 
and self assembly jSS]. The proposed algorithm was modified m to suit 
lattice models like the Lebwohl-Lasher interaction [18], allowing for continu¬ 
ous variation of molecular orientations. It was subsequently augmented with 
the so-called frontier sampling technique [311132] to simulate more complex 
systems like the biaxial medium. The WL sampling is based on effecting a 
convergence of an initial distribution over energy E to the density of states 
(DoS) g{E) of the system iteratively. Frontier sampling technique is an al¬ 
gorithmic guidance, provided in addition to the WL routine, by which the 
system is constrained to visit and sample from low entropic regions. 

The simulations are performed on a cubic lattice of dimensions {L x 
L X L, L = 15and20) with periodic boundary conditions. The biaxial liquid 
crystal molecule on each lattice site interacts with the nearest neighbours 
based on the potential in Eqn.([3]). The uniaxial - biaxial coupling coefficient 
7 on IW is half of the the value on the diagonal IV, for identical A values. 
We denote the arclength of the path OIW as A , given by A = A on segment 
01 , and 

A '= ^(1 + 57 ) 

where 

(1 - 3A) 

on the segment IW. 

The parameters 7 and A were chosen such that we traverse along the path 
IW which amounts to varying the arclength A from 0.33 to 0.747. 
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Table 1: Coordinates of points B1-B7 and Z along the segment IW of the 
essential triangle (Fig.2) 


Point 

7 

A 

A' 

B1 

0.0859 

0.1719 

0.4766 

B2 

0.1405 

0.14658 

0.5666 

B3 

0.1663 

0.1116 

0.6105 

B4 

0.2045 

0.0606 

0.674 

Z 

0.2149 

0.0467 

0.691 

B5 

0.2253 

0.0328 

0.709 

B 6 

0.2440 

0.0079 

0.740 

B7 

0.2482 

0.0024 

0.747 


The Table [T] lists the valnes of ( 7 , A) and corresponding arc lengths A' 
at some designated points for ready reference. Fig. 1(b) shows the location 
of these points schematically inside the essential triangle. Simnlations were 
carried ont for nearly 40 points on IW nsing modified Wang-Landan (WL) 
algorithm angmented by frontier sampling. At each valne, the g(E) which is 
the estimate of the density of states was obtained and an entropic ensemble 
of 10 ^ microstates was generated, by making an effectively nniform random 
walk in energy space gnided by the DoS. Eqnilibrinm ensembles at any de¬ 
sired (rednced) temperatnre (T ) are conseqnently extracted by a snitable 
reweighting procednre [3T1 09] and the average valnes of physical properties 
are calcnlated. The representative free energy F, as a fnnction of the energy 
of the system, as well as of the two dominant order parameters (nniaxial and 
biaxial orders) is compnted from the DoS and the microcanonical energy, - 
both available as a fnnction of bin nnmber in the entropic ensemble. 

The physical parameters of interest in this system, calcnlated at each A , 
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are the average energy < E >, specific heat < Cy >, energy cumulant V 4 
(= 1— < > /(3 < >^)) which is a measure of the kurtosis [SO], the 

four order parameters of the phase calculated according to |23l |5T] and their 
susceptibilities. These are the uniaxial order < i?QQ > (along the primary 
director), the phase biaxiality < i? 2 o >5 molecular contribution to 

the biaxiality of the medium < R 22 >, and the contribution to uniaxial order 
from the molecular minor axes < Rq 2 >■ 

The averages are computed at a temperature resolution of 0.002 units in 
the temperature range [0.05, 2.05]. Statistical errors in different observables 
are estimated over ensembles comprising a minimum of 5 x 10^ microstates, 
and these are compared with several such equilibrium ensembles at the same 
( 7 , A) value, but initiating the random walk from different arbitrary points 
in the configuration space. We find the relative errors in energies are 1 in 
10 ®, while those in the estimation of the order parameters are 1 in 10 ^. 


4 Results and Discussion 

WL simulations were carried out at 40 values of A' on the segment IW, where 
the arc length A ranges from 0.33 to 0.75. Temperature variation of the 
specific heat, and the two order parameters (i?oo R 22 ) in different ranges 
of A , covering the segment IZ (in the attractive region of the interaction 
Hamiltonian), are presented in Figs. [2]-HI The corresponding data along the 
segment ZW (in the repulsive region) is presented in FigJTTl 


4.1 Segment IZ: Range of A = (0.33 - 0.691) 

It is noted from Fig. 2(a) that for all values of A' in the range 0.33 - 0.455, 
a single transition peak is observed in the specific heat profile. As the bi¬ 
axial system is cooled from the high temperature isotropic phase, a direct 
I — Nb transition takes place and the order profiles in Fig. 2(b) reflect the 
behaviour. Fig. 3(a) depicts the splitting of this transition into two transi¬ 
tions for higher values of A in the range 0.462 to 0.610. It is observed that 
lower transition temperature T 2 is progressively depressed with increase in 
A value, as compared to the higher transition temperature peak at Ti. The 
variation of the order profiles in this region, shown in Fig. 3(b), reveals an 
intervening phase which is not strictly uniaxial since the system exhibits a 
low value of R 22 at the onset of the transition. By performing simulations 
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(a) 


(b) 


Figure 2: Comparison of (a) specific heat (b) order profiles for values of A' 
from 0.33 to 0.455 
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Figure 3: Comparison of (a) Specific heat (b) order parameter prohles for 
values of X' from 0.463 to 0.610 






















































at different sizes (L=10, 15, 20), the possibility that this could be a hnite 
size effect is ruled out. (It may be noted that for these system sizes a pure 
uniaxial phase condenses on the A-axis). The notable difference in the case 
of path IW, relative to IV, is that the degree of biaxiality (value of i? 22 ) 
remains fairly independent of temperature, and the degree is the same for all 
subsequent values of A beyond this threshold, until interrupted by a second 
low temperature transition leading to an onset of appreciable biaxial order 

(FigJSD. 




(a) 


(b) 


Figure 4: Comparison of (a) specihc heat (b) order parameter prohles for 
values of A' from 0.639 to 0.691 


The results of simulation for A' values in the range 0.639 - 0.691 are 
depicted in Fig. HI Though these points lie very close to the parabola, they 
are still in the attractive region for the interaction Hamiltonian. It is observed 
from Fig. 4(a) that the second specific heat peak shifts progressively to lower 
temperatures as the value of A increases. Corresponding variations in the 
order parameters, shown in Fig, 4(b) confirm the shift of the second transition 
temperature to lower values. However it is observed that the equilibrium 
averages of order parameters in this region of A are not as smooth, and show 
discernible fluctuations in the low temperature biaxial phase. 

It is interesting to note that the intermediate phase persists to have a 
small degree of biaxial order (~ 0.05) which (a) is not a hnite size effect; 
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(b) is fairly independent of temperatures within the liquid crystal phase; and 

(c) does not depend on the values of A . This phase with temperature de¬ 
pendence typical of uniaxial order, but having a small and constant biaxial 
symmetry (< 0.05)) , is designated as N^i phase in our notation. On sub¬ 
sequent lowering of temperature from this phase, the biaxial order increases 
rapidly at the second transition at T 2 and the lower temperature phase has 
macroscopically observable biaxial order, for all values of A . 



Figure 5: Susceptibilities (y's) of the order parameters for values of A' : 
0.462 - 0.610 

The susceptibility profiles of the order parameter in this region are de¬ 
picted in Fig. m It is observed that the i ?22 susceptibility starts increasing 
in the intermediate phase before showing a peak at the low temperature 
transition at T 2 . 

The fourth order energy cumulant (V4) data obtained along the path OIW 
are shown in Figs. [ 6 ]- El It is observed that the I — Nb transition remains 
strongly hrst order for values of A from 0.345 to 0.45. In the range of A 
from 0.463 to 0.691 ( i.e upto the point Z in FigE]), the high temperature 
transition at Ti from the isotropic phase (I) to the ordered phase shows 
a first order nature. Subsequently, the low temperature Njji — Nb transition 
seems to change gradually from first order to continuous nature, as seen from 
V 4 prohles in Figsj 6 ]andEl 
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Figure 6: energy cumulant V 4 for values of A': 0.345 - 0.45 


An analysis of the above simulation data leads to the proposal of a phase 
diagram along the path IW, shown in FiglO] We could report the L = 20 
data only upto the value A = 0.709, as beyond this value (which falls into 
the partly repulsive region under the parabola) the computational times for 
the convergence of DoS are impractical. We observe from the temperature 
variation of order parameters that the growth of biaxial order appears to 
be progressively inhibited as A value increases within the attractive region, 
and enters the party repulsive region on crossing the parabola at the point Z. 
The free energy prohles, plotted as a function of energy and order parameters 
(computed from the DOS data), reflect the rationale for the impediments for 
the growth of the biaxial order as the base of the triangle OIW is reached. 

The free energy curves obtained for A = 0.610 (Bsin the attractive re¬ 
gion) are shown in Fig. [TOl These curves depict the smooth variation of free 
energy as a function of energy and uniaxial order parameter. However, its 
variation with respect to R 22 shows a small sharp well, (the edge being lo¬ 


cated at R 22 — 0.02), and the family of curves in Fig, 10(c), as a function of 
temperature, shows that it required signihcant variation of temperature be¬ 
fore the system could shift its free energy minimum away from this restricted 
region. It appears that during this temperature range, the system accesses 
microstates with rather small but nonzero degree of biaxiality, constrained 
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Figure 7: Energy cumulant V 4 for values of A' :0.463 - 0.610 
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Figure 8: energy cumulant V 4 for values of A' : 0.64 - 0.691 



A.' 


Figure 9: Phase diagram inside the essential triangle along path IW 
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(c) 

Figure 10: Representative free energy plotted as a function of (a) energy 
(b)i?gQ (c) R 22 at the point B3 (A' = 0.610 ) 
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Figure 11: Representative free energy plotted as a function of (a) energy 
(b)i?oo (c) -R 22 the point B4 (A' = 0.674) 


however by free energy barriers to attain higher degree of biaxiality for con¬ 
siderable range of temperature. This circumstance seems to be manifesting 
as a corresponding curious variation of R 22 at A = 0.610, as in Fig.??. We 
hnd that the development of such free energy barriers (with respect to R 22 at 
low values) and the requirement of the system to cool sufficiently to overcome 
them before accessing higher macroscopically observable values, is generic. 
All the data collected in this region supports and corroborates the simulated 
order parameter prohles reported in the previous figures. It is very interesting 
that such barriers are exhibited only along the path of biaxial order, but not 
along energy or uniaxial order. This implies a complex free energy surface 
in the 2-d space of order parameter, offering initial barriers to a significant 
development of biaxial order, until the system is sufficiently cooled. Figs. [TT] 
-m demonstrate this view point. 
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Figure 12: Representative free energy plotted as a function of (a) energy (b) 
Rqq (c) R 22 point Z (A' = 0.691) 
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Figure 13: Representative free energy plotted as a function of (a) energy (b) 
Rqq (c) R 22 point B5 (A' = 0.709) 
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(a) (b) 

Figure 14: Comparison of (a) specific heat (b) order parameter profiles for 
values of A' in the range 0.709 - 0.747 (L=15) 


4.2 Segment ZW: Range of A = (0.691 - 0.747) 

We now present data obtained beyond the point Z (Fig.??). In this region, 
the biaxial-biaxial tensorial coupling term A ^ 0 asymptotically, leading 
to a special case of the interaction Hamiltonian. The case for A = 0 was 
studied earlier through simulations m It was found that in the absence of 
the biaxial-biaxial interaction term, only a uniaxial phase could be obtained 
on condensation from the isotropic phase. We present here the simulation 
results in the case of A —)■ 0. The mean field analysis predicts that the 
Hamiltonian is partly repulsive in this region and excluded volume effects 
play a major role [8]. Due to the constraints imposed by computational 
time, we could obtain data in this range of A only for a smaller system, 
with L=15. (instead of L=20, as in the earlier case). The specific heat and 
order parameter profiles are depicted in Figs. 14(a) and 14(b) The energy 
cumulants 14 are shown in Fig. [151 It may be observed that the specific 
heat profiles show evidences of two transitions. The order parameter profiles 
depict the onset and growth of uniaxial order at Ti for all values of A . The 
biaxial order parameter increases at T 2 (in the biaxial phase) for A = 0.709, 
but remains close to zero for A = 0.740 and 0.747. This behaviour is as 
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Figure 15: Energy cumulant V 4 for values of A' in the range 0.709 - 0.747 
(L=15) 
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Figure 16: Representative free energy plotted as a function of (a) energy 
(b)i?go (c) RI 2 at point B6 (A' = 0.740) 


expected from mean field considerations at such values of A', very close to 


confirms the presence of barriers for the growth of biaxial order at points close 
to the base OW. The biaxial state is obviously not stable at such parameter 
points of the Hamiltonian. 


the base OW. The free energy plots for A® = 0.740 are shown in Figs. 16(a) 


16(c) The free energy variation with respect to R 22 = 0.740 again 


5 Conclusions 

In this paper, the interior of essential triangle is explored with entropic sam¬ 
pling method along a trajectory IW, a line drawn from apex I of the triangle 
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to the mid point W of the base OV (Fig. [T]). We refer to the arc length OIW 
as A , for purposes of discussion. As per MF treatment this line cuts both the 
trajectories of C 2 C 3 and C 1 C 3 . Thus the phase sequences along this line IW 
should be qualitatively similar to that on the A-axis. In particular the direct 
transition from isotropic to biaxial phase is expected to be interrupted by a 
uniaxial nematic phase beyond a value as 7 increases, and the temperature 
range of the uniaxial phase should progressively increase, suppressing the 
second transition temperature, till the line cuts the parabola, at Z (Fig{T]). 
Results from mean held treatment in the partly repulsive region on this line 
IW are not available for direct comparison, even though it is established that 
biaxial phase would not be stable at the point W PISS]. The phase sequences 
in the present work qualitatively follow this scenario, but with a curious devi¬ 
ation. The intervening ’’uniaxial ” phase Nu is not strictly devoid of biaxial 
symmetry. Indeed all along the line, beyond K (Fig.??), and upto point 
Z, the onset of the uniaxial order is invariably accompanied by a small, but 
unmistakable, development of biaxial symmetry. We thus refer to this phase 
as Njj !, to make this subtle distinction. This small degree of biaxiality of the 
Njj! phase is temperature independent within that phase, and is also fairly 
independent of its location in the trajectory beyond K . An examination of 
the free energy prohles, drawn as a function of both the major order param¬ 
eters, show interesting features: while the free energy curves show smooth 
variation of the minima with respect to Rqq as the temperature is varied, 
the case of R 22 qualitatively different. These prohles exhibit free energy 
barriers at low values of which could be overcome (thereby pushing the 
system to access regions of higher and discernible order), only after these ini¬ 
tial barriers could be overcome on considerable cooling. Thus these results 
show a complex free energy surface that develops with decrease of tempera¬ 
ture on a typical trajectory inside the triangle. It appears that development 
of a Njj! phase with a small biaxial order (< 0.05) is expected, and the degree 
of this symmetry is restricted by the free energy barriers till the system is 
permitted to access these regions of biaxial order. Given that such barriers 
are strongly dependent on the size of the system, it is a plausible conjecture 
to suggest that in real systems these barriers are not readily overcome (or 
equivalently, requires signihcant cooling of the medium), and hence their bi¬ 
axial order appears to be restricted inherently. Under such circumstances 
requiring wider temperature ranges to overcome barriers, real systems may 
have other competing interactions (like translational degrees, influencing the 
phase sequence qualitatively differently , e.g layer formation). Deviations of 
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real systems from MF predictions [S3] could perhaps be understood in these 

terms. 
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